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Abstract: 
In this paper we have suggested two classes of estimators for population median My of the study character 
Y using information on two auxiliary characters X and Z in double sampling. It has been shown that the 
suggested classes of estimators are more efficient than the one suggested by Singh et al (2001). Estimators 
based on estimated optimum values have been also considered with their properties. The optimum values 
of the first phase and second phase sample sizes are also obtained for the fixed cost of survey. 
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1. INTRODUCTION 


In survey sampling, statisticians often come across the study of variables which have highly skewed 
distributions, such as income, expenditure etc. In such situations, the estimation of median deserves special 
attention. Kuk and Mak (1989) are the first to introduce the estimation of population median of the study 
variate Y using auxiliary information in survey sampling. Francisco and Fuller (1991) have also 
considered the problem of estimation of the median as part of the estimation of a finite population 
distribution function. Later Singh et al (2001) have dealt extensively with the problem of estimation of 
median using auxiliary information on an auxiliary variate in two phase sampling. 


Consider a finite population U={1,2,...,i,...,N}. Let Y and X be the variable for study and auxiliary 
variable, taking values Y; and X; respectively for the i-th unit. When the two variables are strongly related 
but no information is available on the population median Mx of X, we seek to estimate the population 
median My of Y from a sample S,,, obtained through a two-phase selection. Permitting simple random 
sampling without replacement (SRSWOR) design in each phase, the two-phase sampling scheme will be as 
follows: 


(i) The first phase sample S,(S,cU) of fixed size n is drawn to observe only X in order to 
furnish an estimate of Mx. 


(ii) Given S,, the second phase sample S,,(S,,CS,) of fixed size m is drawn to observe Y 
only. 


Assuming that the median Mx of the variable X is known, Kuk and Mak (1989) suggested a ratio estimator 
for the population median My of Y as 
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(1.1) 


where M y and M y are the sample estimators of My and Mx respectively based on a sample S,, of size 


m. Suppose that yi), y2), -.-. Yim) are the y values of sample units in ascending order. Further, let t be an 
integer such that Ya) < My <Y¢+1) and let p=t/m be the proportion of Y, values in the sample that are less 
than or equal to the median value My, an unknown population parameter. If p is a predictor of p, the 


sample median M y can be written in terms of quantities as O, ( Pp) where Pp = 0.5. Kuk and Mak 
(1989) define a matrix of proportions (Pjj(x,y)) as 


Y < My Y > My Total 
X < Mx Pilly) Poi(Xy) P(X y) 
x> Mx P1(x,y) P»9(x,y) P.2(x,y) 
Total P,-(x,y) P3-(x,y) 1 
and a position estimator of My given by 
mu,” = 0,(p;) (1.2) 
ae ee Gas y) , (=m) Pale, 2) 
m P(x, y) D»(x,y) 
s of Paley + (m =m) Pa, 2 
m 


with Dp ij (x, ¥) being the sample analogues of the Pj(x,y) obtained from the population and m, the number 


of units in S,, with X < Mx. 


Let F; va(¥) and F vg(¥) denote the proportion of units in the sample S,, with X < Mx, and X>Mx, 


respectively that have Y values less than or equal to y. Then for estimating My, Kuk and Mak (1989) 
suggested the 'stratification estimator' as 


M1, =infly: F, 0.5} ae 
ee kore 
where Fy (y) = a Ba +F,°| 


It is to be noted that the estimators defined in (1.1), (1.2) and (1.3) are based on prior knowledge of the 
median My of the auxiliary character X. In many situations of practical importance the population median 
Mx of X may not be known. This led Singh ef a/ (2001) to discuss the problem of estimating the 
population median My in double sampling and suggested an analogous ratio estimator as 
, ~ M! 
M,, =M,— (1.4) 


X 
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where M . is sample median based on first phase sample S,. 


Sometimes even if Mx is unknown, information on a second auxiliary variable Z, closely related to X but 
compared X remotely related to Y, is available on all units of the population. This type of situation has 
been briefly discussed by, among others, Chand (1975), Kiregyera (1980, 84), Srivenkataramana and Tracy 
(1989), Sahoo and Sahoo (1993) and Singh (1993). Let Mz be the known population median of Z. 
Defining 


e) = ¥ _1|,e,=|—*-l1he, = *__1],e,) —4-1|ande, = 2 
My, M, M, M, M, 


such that E(e,)=0 and | ex | <1 for k=0,1,2,3; where M > and M : are the sample median estimators based 
on second phase sample S,, and first phase sample S,. Let us define the following two new matrices as 


Z<M, Z> Mz Total 
xX < Mx Pi1(x,z) P>1(X,Z) P.1(x,z) 
X > Mx Pi2(x,Z) P39(X,Z) P.2(X,2Z) 

Total P,-(x,z) P,-(x,z) 1 

and 

Z<Mz Z>Mz Total 
Y <My Piily,Z) Poily,2Z) P.i(y.Z) 
Y >My Pialy,zZ) P20 y,Z) P.o(y,Z) 

Total P,-(y,z) P2-(y,zZ) 1 


Using results given in the Appendix-1, to the first order of approximation, we have 


E (ex?) = ("3") amy" Myfy(M)}, 
B¢e?)= (SJ) amy" MxM}, 
(er) = (YP) any MxM}, 
(ex) = (53°) 4m)" MetM}° 
(ex?) = ("Y) ny (Met) }”, 


E(eoe}) = Mn) (4m)" {4P 11(x,y)- 1} {MxMyfx(Mx)fy(My)} "1, 


M2) (4m {4P11(y)-1} {(MxMyfic(Mso fy(My)}" 
B(ee:) = (SIP) (amy (4P fy.2)-1} {MyMrf(My (M2) 1 


E(eeq) = xn) (4n)! {4P11(y,z)- 1} {MyMefy(My)f2(Mz)}", 


B(e.e:) = (92) ny MyM}, 


E(e\e3) = Mn) (4m) ' {4P11(x,z)-1} {MxMzfx(Mx)fz(Mz)}", 
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<=} (4ny! {4P13(x,z)-1} {MxMzfx(Mx)f,(Mz)} 
S) (An) {4P 11(x,2)-1} {MxMrfx(Mx)£(Mz)} 


E(e.€4) = a (4n)'{4P11(x,z)-1} {MxMzfx(Mx)fz(Mz)}", 


E(e3ea) = (= } (4n)(f(Mz)Mz)” 


where it is assumed that as No the distribution of the trivariate variable (X,Y,Z) approaches a continuous 
distribution with marginal densities fx(x), fy(y) and fz(z) for X, Y and Z respectively. This assumption 
holds in particular under a superpopulation model framework, treating the values of (X, Y, Z) in the 
population as a realization of N independent observations from a continuous distribution. We also assume 
that fy(My), fx(Mx) and f7(Mz) are positive. 


Under these conditions, the sample median M , is consistent and asymptotically normal (Gross, 1980) with 
mean My and variance 


In this paper we have suggested a class of estimators for My using information on two auxiliary variables X 
and Z in double sampling and analyzes its properties. 


2. SUGGESTED CLASS OF ESTIMATORS 
Motivated by Srivastava (1971), we suggest a class of estimators of My of Y as 


g = (1, a7, =M,e(u,v)} (2.1) 


A ~ 


x v=—~ and g(u,v) is a function of u and v such that g(1,1)=1 and such that it satisfies 


where uv = ee 


Xx Zz 
the following conditions. 


1. Whatever be the samples (S,, and S,,) chosen, let (u,v) assume values in a closed convex sub- 
space, P, of the two dimensional real space containing the point (1,1). 


2. The function g(u,v) is continuous in P, such that g(1,1)=1. 
3. The first and second order partial derivatives of g(u,v) exist and are also continuous in P. 


Expanding g(u,v) about the point (1,1) in a second order Taylor's series and taking expectations, it is found 
that 


E(7,®)= My +00) 
so the bias is of order n7?. 


Using a first order Taylor's series expansion around the point (1,1) and noting that g(1,1)=1, we have 
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mM, = M, [1 + €o =F (e, =: e, gy (1,1)+ C482 (1,1) + O(n) 
or 
(u (2) —M,)= M, le, + (e, —e, )g, (11) + €482 (1) 22) 


where g)(1,1) and g,(1,1) denote first order partial derivatives of g(u,v) with respect to u and v respectively 
around the point (1,1). 


Squaring both sides in (2.2) and then taking expectations, we get the variance of M oe 


rar," (LE }+(2-2)ar(2-2 ha), ea 


A= Baan Z, on g,(11)+ 2(4P,,(x,¥)— i (2.4) 


to the first degree 
of approximation, as 


where 


B= Earnie es (| (ee) gy (11) + 2(4P,,(y,z)—- i (2.5) 


(g) 


The variance of M y in (2.3) is minimized for 


g(1) = oe a ita an (oy)-1 
(M 


ext) =| $2 


(2.6) 


Thus the resulting (minimum) variance of M ys) is given by 


min. Vary") = (; Ee 7 AP bsa) I (Fa J4Aa0.2)-9] 


Now, we proved the following theorem. 


Theorem 2.1 - Up to terms of order n"', 


vari) era) Gefen) (Fy Mabe) a | 


with equality holding if 
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g (1) = {eet Nantes) -1) 
g2(ll) = {ee an, (y,2)-1) 


(g) 


It is interesting to note that the lower bound of the variance of MM . ®” at (2.1) is the variance of the linear 


regression estimator 


n 


M1, =M, +4,(My -M,)+d,(M, ~M') (2.8) 


with By (x,y) and By ly,z) being the sample analogues of the Di (x,y) and Pi (9,2) respectively 
and f y (7 . } Z 7 (mM ze ) and 4 ; (mM 5) can be obtained by following Silverman (1986). 
Any parametric function g(u,v) satisfying the conditions (1), (2) and (3) can generate an asymptotically 


acceptable estimator. The class of such estimators are large. The following simple functions g(u,v) give 
even estimators of the class 


i ap (2 _ 1+a(uv-1) 
Ou,v)=u ‘4 ee (u, = Gai)’ 
g(u,v)=1+a(w-1)+ B(-D,g (wv) = {l-a(w-1)- BI 
g(u,v) = wu? +w,v?,w, tw, =1 

g )=au+(l-a)v? ,g(u,v) = exp{a(w—1)+ B(v-1)} 


Let the seven estimators generated by g")(u,v) be denoted by M a? =M : gu, v), (i =1 to7). It is 


easily seen that the optimum values of the parameters o,8,w;(i-1,2) are given by the right hand sides of 
(2.6). 


3. A WIDER CLASS OF ESTIMATORS 


The class of estimators (2.1) does not include the estimator 
M yy =M, +d, (71 -M,)+d,(M, -M)(d,,4,) 


being constants. 
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However, it is easily shown that if we consider a class of estimators wider than (2.1), defined by 
M © _ a(x u,v) (3.1) 
Y el yous : 


of My, where G(-) is a function of M, , u and v such that G(M, Ail) = M, and G, (M, 1,1) =], 
G, (M e 1,1) denoting the first partial derivative of G(-) with respect to M if 


Proceeding as in Section 2 it is easily seen that the bias of M Mo) is of the order n' and up to this order of 


(@) 


terms, the variance of M, ’ is given by 


varlits”)- aye eae AACA) 
G,(M, DO Reatale (M,,,1,1)+ 2(4P, (x,y) 7 


fi 1) Al) If fy) 
cea cacal 


n N f,(M,)M, M,f,(M, 


flatts at4R 42)-0h 
(3.2) 


where Go(My1,1) and G3(My1,1) denote the first partial derivatives of u and v respectively around the point 
(My,(1,1). 


: me (cs NV a 
The variance of MZ ut ) is minimized for 


6,(a, aa)=-{ Mate) Var, (s.9)-1 


G, (uM, 1,1) = {Meee lan (y, z) = 1) 


(3.3) 


Substitution of (3.3) in (3.2) yields the minimum variance of M Je) as 


in. Vanity) = mal ]-( 2 fara lsn) th (2 Jarre) 
= min Var," 
(3.4) 


Thus we established the following theorem. Theorem 3.1 - Up to terms of order n", 
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var(ity 2 (4-2) (dar (s.n)-0) (4-2 Jari (2)-1) | 


a(f,(My)) Lm NJ \im n 


with equality holding if 


6,(M, td) =| A Yar, (sy) 


G, (uM, 1,1) a {Ata Nan, (y, z) a 1) 


: : me i . . ~ (G 
If the information on second auxiliary variable z is not used, then the class of estimators MZ a reduces to 
the class of estimators of My as 


M7, = H(M,,u) (3.5) 
where H(i, u) is a function of (x7, u) such that H(M, 1) = M, and H, (M, 1) =], 
i (M, 1) = oH(') . The estimator mM, is reported by Singh ef al (2001). 


¥ 4(My,1) 


jo : “ (H spate Atte thts oct 
The minimum variance of M f ) to the first degree of approximation is given by 


min Vary" }=@ (4 2)-(2 2 Var,Gxy)-1) 66 


af (My Lon) in 


From (3.4) and (3.6) we have 


1 \) 1 
n N)Af,(M,)) 


minVar(Mt,” \s min. Var(My°) ) = 5 (4P,(y,z)—-1)° G7) 


ae 3 : ~ (G). ‘ : 
which is always positive. Thus the proposed class of estimators MM a ) is more efficient than the estimator 


M,” considered by Singh et al (2001). 


4. ESTIMATOR BASED ON ESTIMATED OPTIMUM VALUES 


We denote 
_Mxfy (Mx) 4P =] 
“a M,fAM,) | ney bi 
_MifWa\ap (5)- | 
Qs, = ParaVTE Gale ) 1) 
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In practice the optimum values of g)(1,1)(=-a,) and go(1,1)(=-a2) are not known. Then we use to find out 
their sample estimates from the data at hand. Estimators of optimum value of g)(1,1) and g,(1,1) are given 
as 


(4.2) 


where 


(4.3) 


Now following the procedure discussed in Singh and Singh (19xx) and Srivastava and Jhajj (1983), we 
define the following class of estimators of My (based on estimated optimum) as 


mM, =M,g*(u,v,d,,0,) (4.4) 


where g*(-) is a function of (u,v, ,,@,) such that 
as (1,1,0.,0., ) = 


sitieo,)= 220) =n 

g2(La,.0%)= 1b ax = 

g:(1La,,0,)= am & 0 
11,0, 0, 

g.(L1o,,0,)= 7 | =0 
ade 


and such that it satisfies the following conditions: 


1. Whatever be the samples (S, and S,,) chosen, let u, V,0 4 assume values in a closed convex sub- 
space, S, of the four dimensional real space containing the point (1,1,01,02). 


2. The function g*(u,v, 0, 2) continuous in S. 


3. The first and second order partial derivatives of g * (u, V0, Kol 5 ) exst. and are also continuous in 
S. 


Under the above conditions, it can be shown that 


E(m,")=M, +0(n*) 
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and to the first degree of approximation, the variance of M yf? is given by 
Var(a7, ) = min. Var(My* ] 


where min. Var(47,“’) is given in (2.7). 


A wider class of estimators of My based on estimated optimum values is defined by 


2 (G*) a s Axe AK 
M,' =G*\M,,u,v,a, 0, 
where 


ap Mahe) ias 6,y)-1) 


a3 = Se Hi Mp0) 


are the estimates of 


and G*(-) is a function of (7, ,U,V,0, a3) such that 


G*(M, 11a; a3)=M, 


G?(M, 11,a*,03)= se () 
Y 


. = JOG XC 
G3(M, 1la% a5 a6") 


=1 


(1, JLLa* a3) 


=O, 


(My A.Laf 03) 
aG *(-) 


ov (uy 1,1,07a3) 


G3(M ,1,La°,05)= =O; 


G;(M 1,10; 03 )=— =0 


(uy JLo; a3) 


(4.5) 


(4.6) 


(4.7) 


(4.8) 


Under these conditions it can be easily shown that 


n 


E(m,'°)=M, +00) 
and to the first degree of approximation, the variance of M at is given by 

Var a ) = min. Var(Z ) (4.9) 
where min. Var(AZ =") is given in (3.4). 


a * 
It is to be mentioned that a large number of estimators can be generated from the classes M eo ) and 


“ (G 
M based on estimated optimum values. 


5. EFFICIENCY OF THE SUGGESTED CLASS OF ESTIMATORS FOR FIXED COST 


The appropriate estimator based on on single-phase sampling without using any auxiliary variable is M, , 
whose variance is given by 


n 1 1 1 
vali) Saga 2 


In case when we do not use any auxiliary character then the cost function is of the form Co-mC,, where Co 
and C, are total cost and cost per unit of collecting information on the character Y. 


The optimum value of the variance for the fixed cost Co is given by 


ont va fe ) =i [2 - *)| (52) 
where 
V, nt oe (5.3) 
4(f,(M,)) 
When we use one auxiliary character X then the cost function is given by 
Cy, =Gm+Cy,n, (5.4) 


where C; is the cost per unit of collecting information on the auxiliary character Z. 


(H) 


The optimum sample sizes under (5.4) for which the minimum variance of M °° is optimum, are 
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Cy (V, -V,)IC, 
Mca -V,)c, +WV.C | 


(5.5) 


11 


Coaivia Cs 
WW, -V,)c, + IVC] 


opt 


where V;=V0(4P);(x,y)-1)”. 


: ; : oe ; : ~ (H). 
Putting these optimum values of m and n in the minimum variance expression of M A ) in (3.6), we get 


the optimum min. Var(\7 Ae) ) as 


2 
nen the (Vc, +\V.C,) v 
Opt [min.Var(i7 ut) J = es ! as : (5.7) 
C, N 
Similarly, when we use an additional character Z then the cost function is given by 
C, =C,m+(C, +C,)n (5.8) 
where C; is the cost per unit of collecting information on character Z. 
It is assumed that C;>C,>C3. The optimum values of m and n for fixed cost Cp which minimizes the 
minimum variance of M y®) (ora 7) (2.7) (or (3.4)) are given by 
Co, -“Y/C, (5.9) 
opt — : 
NW, =H, +I +O, -7)] 
Col, -V, iC, +C, 
ee (5.10) 


WO VIC, +(C, +, -%)| 
where V>=V(4P 1:(y,z)-1)”. 


The optimum variance of M a (ori oe ) corresponding to optimal two-phase sampling strategy is 


= a 2 
Opt|min. Var(a7 1 or min.Var(7 } = IMs . ie s ug = M Z ) ~ 
0 
(5.11) 


Assuming large N, the proposed two phase sampling strategy would be profitable over single phase 
sampling so long as 


lOpt. Var(i7 y ] > Opt|min.Var(17 ve ) Jor min. Var(\7 wo J 
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Sacer: < Wo = WVo=V, 
C, VM -V2 


When N is large, the proposed two phase sampling is more efficient than that Singh et a/ (2001) strategy if 


(5.12) 


Opt[min. Varli7 Me br min. Var(/ A J i Opt[min. Var(7 mal) 


+ 
ae ee 


i. (5.13) 
C, V, i V, 
6. GENERALIZED CLASS OF ESTIMATORS 
We suggest a class of estimators of My as 
~ oF). yy (F y 
3=(M, a, = F(A, ,u,v, w}} (6.1) 


where u =M, /M'., Vv =M', /M,, w=M, /M , and the function F(-) assumes a value in a 
bounded closed convex subset WCR,, which contains the point (My,1,1,1)=T and is such that 


F(T)=My=>F |(T)=1, F,(T) denoting the first order partial derivative of F(-) with respect to M y around the 
point T=(My,1,1,1). Using a first order Taylor's series expansion around the point T, we get 


M, = F(T) +(M, =M, )R(T)+ U-DE(T) + (v- DET) + w-DE,(1) + 000") 
(6.2) 


where F(T), F3(T) and F,(T) denote the first order partial derivatives of F’ (7 yoU,V, wv] with respect to u, 


v and w around the point T respectively. Under the assumption that F(T)=My and F,(T)=1, we have the 
following theorem. 


Theorem 6.1. Any estimator in 3 is asymptotically unbiased and normal. 


Proof: Following Kuk and Mak (1989), let Py, Px and Pz denote the proportion of Y, X and Z values 
respectively for which Y<My, X<Mx and Z<M7; then we have 


My My = ap tl-2F,)+0,(" 
We My =z yl-2F)+0,(2") 
Me a5 Ge je tee") 
= = Fpuayt tel ele) 
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M,-M, sary 2P)+ (9) 


Using these expressions in (6.2), we get the required results. 


Expression (6.2) can be rewritten as 


M7, -M, =(M, -M,)+(u-DF (7) + (v-DF(T)+ (w-DF,(T) 


mM, —M, = Mye, +(e, -e, )F,(T) + e,F,(T) +e; F, (7) (6.3) 


Squaring both sides of (6.3) and then taking expectation, we get the variance of M a) to the first degree 
of approximation, as 


vali le ap ( ata i Cow} = 


4-|to( fed) ear satin O22) i Alte) len 


M,f,(M, 


{itt em auncen-nee 
A -( f(y) [wr M,) a )+ ( 11 (X,Y) ) a ) 
Myf: (My) + 2(4P,, (x, Z) (a2 jJenvesn 


f,(My) 


M,f,(M, A (T)+ 2(4P,,(y,Z) —1)F, (T) 


fr(My) 
GOO 


The Var i) at (6.4) is minimized for 
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[(4P (x,y) -(4F i (.2)- 4A (y2)-U Mate (Mx) 


F,(T) = ; 
= (1 (4P,,(x,2)—1) fAMy) 
= —a, (say) 
(6.5) 
rer) = Paz) MGR (sy) =) (4 (v2) AR (2) Ma fe (M2), 
; [1-(4P,,(x,z)-1)"] f, (My) 
= —a, (say) 
yr) = —WaPaly.z)=)= GP, (e.y)=D4Pa(s2)-D Mz fe(M 2), 
, [1-(4P,(x,z)-1)"] f,(M,) 
= —a, (say) 


Thus the resulting (minimum) variance of M se) is given by 


minVar(7 a ) = 


1 | D 
m n)4(f,(M,)) 1-[4P,(x.z)-1| 


(6.6) 


where 
D =(|(4P, (y,z)-1)-(47, (x, »)-1(4F, (xz)-)] (6.7) 
and min. Var(i7,] is given in (3.4) 


Expression (6.6) clearly indicates that the proposed class of estimators M a is more efficient than the 


class of estimator M Vo or (7 )] and hence the class of estimators M i suggested by Singh et al 
(2001) and the estimator M y at its optimum conditions. 
The estimator based on estimated optimum values is defined by 

pe =) -M," =F *(M,,u,v,w,4,,4,,4, } (6.8) 


where 
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n 


,, — [4iule.r)-1)- Gauls2)-4u0.2)- DMF.) 


| (iu (2)-1)) f(y) 
a, = (4,,(x,2)— I(4A,, (x »)-1)- (48,,(y.2)- (4A, (%2)- 0] M,f,(M,) 
: 1—(49,,(x,z)-1) f(y) 
_ [(46,,(y.z)-1)—(4P.(x, y)-1) 4 (xs2)-)) Me Fl) 
: 1- (48,,(x,2)-1)"| f,(M,) 


(6.9) 


are the sample estimates of a), a, and a3 given in (6.5) respectively, F*(-) is a function of 


M,,,U,V,W,@,,4y,4, } such that 
*(T#) — 
F*(T*)=My 


= Fer) =F *0) =1 


y (pe 
rr") __, 
Cu. |e 
F,e(r*) =U) cond 
Ov |r 
*/(. 
RAT) =S 0 =a, 


—” 
II 
So 


*(. 
rer) = OE 


1 T* 


F, *(T*) = oF *(:) =.() 


An 


2 T* 


r+) = 9 
a, 


where T* = (My,1,1,1,a1,42,43) 


Under these conditions it can easily be shown that 


E(,")=M, +007) 
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(F*) 


and to the first degree of approximation, the variance of M y is given by 
Var Bae ) = min. Var(\7 oa (6.10) 
where min. Var(\7 ww) ) is given in (6.6). 


Under the cost function (5.8), the optimum values of m and n which minimizes the minimum variance of 


mM,” is (6.6) are given by 


its CyY -V, -V;)/C, 6.11 
iY -V, -W)C, + (U, -% VNC, +3) ie 


CoV -V, =ValC; 
ia LW, a -V,)C, + (V, -V, +V\(C, +C;)] 


where 
D°V, 
V~, = ° (6.12) 
 [@R.G2)-1)] 
for large N, the optimum value of min. Var(\7 ) is given by 
Opt [min.Var(7,”) _ LW, = V, = V; Ie, a VV, = V; + V, Nc, + C; )| (6.13) 


Co 


The proposed two-phase sampling strategy would be profitable over single phase-sampling so long as 


Opt |Var(i7 ; ] > Opt|min.Var(M1,”’ J 


2 
C,+C, |e none 


(6.14) 
m1 Vv V, = V; + V; 
It follows from (5.7) and (6.13) that 
Opt min. Var(i7 eu J < Opt|min. Var(i7 J 
if cI iar et ie Speco a ©, (6.15) 
Vi —V,+V; C, (V, —V, +V ,)C, C, 


for large N. 


Further we note from (5.11) and (6.13) that 
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Opt |min.Var(17 ) ] < Opt|min.Var(17 vf or Me J 


G+, | We—W)-WO—)) 


if 
C, VY. -V +V)-VWV HV; 


(6.16) 
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